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Abstract 

Decoherence effects on quantum and classical dynamics in reactive scattering are examined using 
a Caldeira-Leggett type model. Through a study of dynamics of the collinear H + H2 reaction and 
the transmission over simple one-dimensional barrier potentials, we show that decoherence leads 
to improved agreement between quantum and classical reaction and transmission probabilities, 
primarily by increasing the energy dispersion in a well defined way. Increased potential nonlinearity 
is seen to require larger decoherence in order to attain comparable quantum-classical agreement. 
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I. INTRODUCTION 



Isolated molecular collisions have been the subject of theoretical study over the past 
half century. Amongst the successful tools in the field are classical trajectory approaches, 
whose justification often relies upon the relatively small deBroglie wavelengths in a chemical 
reaction. However, most chemical reactions occur in the condensed phase. Such an environ- 
ment, due to its randomness, tends to decohere flj a quantum system, driving it towards the 
classical limit. Hence, we would anticipate that the utility of classical methods would gain 
further validity from this loss of phase information. In this paper we investigate the extent to 
which quantum effects are lost, in chemical reactions, due to the environment. Specifically, 
we examine, using a Caldeira-Leggett type of approach applied to the paradigmatic 

collinear H + H2 system 5] , the effects of the environment on the two predominant quantum 
effects in reactive scattering: tunneling and resonances. We find, for reasonable values of 
the decoherence, that the quantum dynamics of the (system + environment) approaches 
that of the classical dynamics of the (system + environment). A detailed analysis of the 
tunneling and resonance regimes is also provided through a study of simple one dimensional 
models. Improved decoherence-induced quantum-classical agreement in the reactive cross 
section is found, in both the tunneling and resonance regions, to stem from increasing energy 
dispersion with increasing decoherence strength and increasing interaction time. 

The general issue of the transition from quantum to classical dynamics remains a subject 
of considerable interest and activity. Traditionally, classicality is argued to emerge as h goes 
to zero or equivalently as the system mass goes to infinity. A comprehensive Liouville-based 
approach [6] provides a more complete description of this limiting process for both integrable 
and chaotic systems. Alternatively, decoherence, i.e. , the loss of quantum coherence as 
induced by the interaction of a system with its environment, has been proposed to be 
an essential ingredient for the utility of classical mechanics in the real world [1]. Early 
theoretical studies on model systems successfully demonstrated that decoherence is effective 
in bringing about classical behavior in the macroscopic limit Q|. In addition, our earlier 
study[8] extended this work to Hamiltonian systems, showing the emergence of classical 
behavior due to decoherence in a bound Hamiltonian system. 

De splt e _ papers, and son* appHcatioas a, « UUle - has been 

done on the role of decoherence in inducing classicality in scattering problems, such as 
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chemical reactions. We do so below, via a simple master equation for the reduced density 
matrix p, i.e. the density matrix for the scattering system after averaging over the external 
bath. The essential observable under consideration is the reactive scattering cross section, 
which displays both quantum tunneling and resonance scattering. 



II. DECOHERENCE: FORMULATION AND COMPUTATION 



Consider a scattering event that occurs in the presence of an external environment, focus- 
ing attention on the short time regime wherein decoherence, but not energy transfer between 
the system and environment, occurs. A master equation for p for such a setup has been 
obtained by Joos and Zeh [2] for an isotropic environment that elastically scatters from the 
system of interest with negligible momentum transfer. With the additional assumption that 
the (x\p\x') of interest are those with k\x — x'\ << 1, where k is the typical wavenumber of 
the environment and x is the system coordinate, they obtained an equation for p that is the 

n n 

same as the oft-cited model of Caldeira and Leggett |3J and that of Unruh and Zurek 
when dissipation is ignored. The latter corresponds to ignoring the transfer of momentum 
to the environment, a reasonable assumption for the short time scales associated with the 
collision process. 

For a system with two degrees of freedom, the resultant quantum master equation is 



given as 
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with the system Hamiltonian H and the decoherence constant D. To compare to the classical 

n 

result it is convenient to consider this equation in the Wigner representation |3J: 

§- t P w = Up w + Kp w + Ldp w ; (2) 

Lci = {H, } PB , (3) 

L ° ~ D \dp 2 x + dpi) ' ^ 

where {-}pb denotes the classical Poisson bracket and where p w = p w (p x ,P y ,x,y) denotes 
the Wigner transform of p. Here, The L c i generates the classical dynamics of the closed 



system, and L q generates the quantum corrections to this dynamics, whereas I_d induces 
decoherence in the open system. When D = Eq. (J2J) becomes the quantum Liouville 
equation for the closed system. 

Equation (J2J) is solved using the quantum state diffusion (QSD) approach j^, Q] wherein 
one solves an associated stochastic Schrodinger equation for the state vector given by 

|#> = (-Lh-Y, [Lm- (L m ) M } 2 ) \^)dt + Y[L m -(L m ) m ] \4>)dW m . (6) 

Here the self-adjoint operators L m represent the coupling between the system and environ- 
ment, (L m )\^) = (%l)\L m \ip) I (if}\if}) , and W m is a complex Wiener process jl^. Our specific 
implementation of the QSD method takes the operators L m as L\ = ^-x and L2 = 
Computationally, the unitary component of Eq. © is integrated by the fast Fourier trans- 
form method plus the split-operator ^| , as in a closed quantum system. The terms involv- 
ing L m are integrated using a second order scheme Q]. The dynamics resulting from the 
Stochastic Schrodinger equation (|BJ) . when averaged over many realizations of the Wiener 
process, provides the solution to Eq. (J2J Q| . Hence, the quantum reaction probabilities for 
the open system are obtained by averaging over reaction probabilities for each wavefunction. 

A comparison of these results with the classical result means comparing the quantum 
result of the (system + environment) with the classical result for the (system + environment). 
The latter is obtained by setting L q = in Eq. (j2J, leading to the classical Fokker-Planck 
equation for the classical density p cl = p cl (p x ,p y ,x,y): 

|p cl = L clP cl + L D p d . (7) 

The resultant dynamics is equivalent to the Langevin-Ito equations with Gaussian white 
noise for each degree of freedom, which are solved by combining Monte Carlo sampling from 
the initial distribution and fourth-order Runge-Kutta integration. Results in the absence 
of decoherence (D = 0) were obtained directly by solving Hamilton's equations for the 
classical trajectories. The classical reaction probabilities for both the closed and open cases 
are obtained by counting the number of trajectories that are reactive. 
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III. H + H 2 



A. Computational Features 

We present results for collinear H + H 2 . Isotopic variants of this system were also studied 



22j |. with similar results to those reported here. The Hamiltonian for the collinear H + H2 
system is H = j^(pl + P^) + V(x,y), where x and y are the mass-scaled Jacobi coordinates 
of R (the distance of H from the center of mass of H 2 ) and r (the H 2 internuclear distance), 
respectively, and p x and p y are the corresponding momentum operators. The reduced mass 
H = y[m 1 m 2 m3/ (mi + m 2 + m 3 )] where m.j is the mass of H. The Liu-Siegbahn-Truhlar- 
Horowitz potential energy surface lf| was used for V(x,y). 

The initial wavefunction, placed far from the interaction region, can be written as the 
product of a minimum uncertainty Gaussian wave packet F(x) describing the relative ini- 
tial translational motion times the vibrational eigenfunction <p v {y) for H 2 with vibrational 
quantum number v. That is, ip{x,y,t = 0) = F(x)4> v (y), with 

W \ (n 2N-1/4 r i.x-x f ip x0 x 

F{x) = {2u-i) exp[ — — ]. (8) 

Here 7 is the width of F(x), and Xq and p X Q describe the locations of the initial wavepacket 
in position and momentum, respectively. A large value of 7 results in a spatially delocalized 
wavepacket that is localized in momentum space. The initial total energy, E, is a sum of 
the average kinetic energy of the translational wavepacket and the vibrational energy. In 
the study below we consider the ground vibrational state (v = 0) and vary the translational 
energy. Note that the energy width 5E = J (H 2 ) — (H) 2 of the initial wavepacket, of 
considerable interest later below, is determined by the translational component. Thus, by 
choosing 7 appropriately the wavepacket can be focused on one average energy component 
of interest, E, with a narrow energy width. Specifically, since the initial wavepacket is in 



the asymptotic region [x — > 00) and p x0 is chosen as p x0 = —y2fi(E — E v ), 5E is given by 
5E = y / ' 32^2 + h av^f^ ' w hh vibrational energy E v . Since the smaller the 7, the faster 
the numerical computation, 7 is taken as small as possible, limited by the need to obtain an 
adequate SE that depends on the resolution of energy region studied. Our choice of 7=5.9 
a.u. gives a 5E w 0.068 eV for the initial wavepacket, which is somewhat larger, due to 
computational limitations, than that expected analytically. 
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To examine quantum-classical correspondence requires an appropriate comparison with 
an initial classical ensemble corresponding to the above initial quantum state. For the 
translational component, this is given by the Wigner function corresponding to F(x) in Eq. 
(JSJ). However, there are several possible choices for the vibrational component, w 
used to make the initial quantum and classical conditions as similar as possible 
For example, since the initial state is a vibrational eigenfunction, the corresponding Wigner 
function does not vanish at the classical vibrational turning points, so that some phase points 
have initial classical energies that do not equal the quantized value. The resultant SE for 
the vibrational motion differ for the initial quantum and classical density functions; the 
vibrational energy width of the initial quantum state is zero but that of the classical density 
is hu/2. For the case of H 2 , this equals 0.27 eV. This width is much larger than that of 
the quantum counterpart whose entire energy width, 0.068 eV, comes from the translational 
part. This leads to classical reaction probabilities that are insensitive to decoherence. 

Our results show that for studies of reaction probabilities versus energy, it is crucial to 
have the same initial SE for both the quantum and classical systems. Hence, each classical 
trajectory is chosen to have the fixed quantum vibrational energy, guaranteeing the same 
SE as that of the quantum system. We assume that Ho is described as a Morse oscillator 
and obtain the initial vibrational variables as in Ref. HQ- 

Below, two different values of D are examined. The smaller of the two, D = 2.47 x 10~ 35 
kg-J/s, corresponds to 10~ 2 to 10 _1 of typical D values in solution (a value estimated in 



Appendix). The larger D is 20 times the size of the smaller. Both these values of D are in 



the range extracted for D from a pseudo-realistic model problem 21]. For the open quantum 
system we averaged over 72 realizations for the smaller decoherence and over 288 realizations 
for the larger one. For both the closed and open classical systems 10 4 trajectories are used. 
Below, the notation QM and CM denote the quantum mechanical and classical probabilities, 
respectively, in the absence of decoherence. The symbols QMD and CMD denote results in 
the presence of decoherence. 

We focus below on the effect of decoherence on the total reaction probability P R (E), 
which can be written as 

P R (E) = / J R (t, E)dt = Aty J R (nAt, E), (9) 
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with a chosen time increment At. Here the reaction flux, J R (t, E) is given as: 

J R {t, E) = (h/fM)lm[(ip E (x, V = 3ft, t)\frj> E (x, y = y h t)/dy)), (10) 

where iPe{x, y, t) is the system wavefunction at time t that evolves from an initial wavepacket 
i/jE(x,y,t = 0) that is very narrowly focused on E, and the brackets denote the integration 
over x, obtained by adding all the discrete contributions over the entire range of x along a 
dividing line y = yi in the product region. 

The simulation uses absorbing boundary conditions that damp the wave packet with a 
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sine masking function near the edge of the grid in both reactant and product channels 
We used the following parameters for QM and QMD; the grid starts at (x,y) = (x m in, J/min) 
= (0.90,0.38), the mesh size (Ax,Ay)= (0.15,0.12), the initial position (x , 2/o) = (30, 1.35), 
and yi = 3.02 in atomic units. For QM, At= 0.14 fs, and for QMD At is 0.073 fs for the 
smaller D and 0.012 fs for the larger D. The total time steps are chosen appropriately, 
keeping tf = 287 fs constant. The classical reaction probability is evaluated by counting the 
number of trajectories that cross the dividing line y = yf M = 5 a.u. before the final time 
ty M =290 fs. Numerical checks of quantum and classical calculations simulations at different 
temporal and spatial resolutions were carried out, and results for the closed systems were 



checked against previously computations! 



0,Q,Q. 



B. Computational Results 

Figures Q and El display the classical and quantum reaction probabilities as a function of 
E for the H + H 2 exchange reaction in the closed system and open systems for both values 
of D. For the closed system the results agree well with those of the previous studies Q] and 
one sees large deviations between the classical and quantum probabilities. In particular, note 
the strong quantum resonance dips at E ~ 0.9 eV E ~ 1.2 eV as well as the tunneling that 
are all absent in the classical result. Upon introducing the smaller of the decoherence, the 
resonances are damped out, as seen in Fig. ^ leading to better agreement between quantum 
and classical reaction probabilities. Also notable is that the addition of the decoherence 
in Fig. enhances the reaction probabilities at E < 0.55 eV and suppresses them at 0.55 
eV < E < 0.7 eV, for both the quantum and classical cases. This behavior is even more 
pronounced in the case of the larger D, shown in Fig. |21 and is analyzed in Section lV Bl where 
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tunneling in one dimension is examined in detail. Note from Fig. [I] that the CMD and QMD 
results are still in considerable disagreement. This disagreement disappears with the higher 
D [Fig. |2] where the classical and the quantum results are in very good agreement over the 
entire energy region. However, in this case most of the scattering structure is suppressed: 
the reaction probabilities are elevated at energies below the threshold and reduced above 
the threshold energy, and the Pr curve is now only slightly convex about Pr = 0.5. 

Note that we also examined the case where the Wigner function is used for the vibrational 
component of the initial classical ensemble. In this case the classical reaction probabilities 
showed, in the threshold region, a slower rise with energy than the quantum counterparts 
in both the closed and open systems. This behavior is due to the larger 5E of the initial 
classical density. In this case the resultant insensitivity of the classical reaction probabilities 
to decoherence leads to a persistent discrepancy between quantum and classical reaction 
probabilities, especially in the threshold region. Hence, in this case, basing the classical 
density on the Wigner function provides misleading results. 

Having demonstrated that classical dynamics emerges in the presence of decoherence for 
the H + H2 reaction, it remains necessary to identify the specific mechanism by which this 
occurs in the Caldeira-Leggett type model. To do so, we first show below that the system 
energy dispersion 5E = [(H 2 ) — (H) 2 ] 1 ^ 2 increases in these models as y/D/m. Tunneling 
and resonance behavior for a simpler one dimensional model are then examined to establish 
the role of this increasing SE as the source of the decoherence-induced quantum-classical 
correspondence. We then return to remark on the case of reactive H+H 2 . 



IV. TIME EVOLUTION OF THE ENERGY DISPERSION 5E 

The friction free version of the Caldeira-Leggett master equation [Eq. is known to 
display an increase in system energy due to the D dependent decoherence term (the so-called 
environmental localization term). This energy increase would be compensated for, at far 
longer times, by the transfer of energy back to the environment [12]. More significantly, for 
our purposes, is the fact, shown below, that over similar timescales the square of the energy 
dispersion {SE) 2 = (H 2 ) — (H) 2 increases as D/m 2 , where m is the system mass. To see 



this, we extend the method of Ref. 



12j to the case of energy dispersion. 
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For simplicity, consider a one dimensional case, in which Eq. (0) reduces to 

f t =l[H, P ]-^[ X ,[ X ,p]]. 

For any observable O that does not have explicit time dependence: 

d(O) = l Tr(0p) 
dt dt 



±Tv(0[H,p}) -^Tr[o[x,[x,p}}). 



We note two useful relations: the cyclic invariance of the trace, which gives 

Tr(A[B,C]) = Tr([A,B]C), 



(11) 



(12) 



(13) 



and the identity 



[AJ(B)] = df/dB[A,B], 



(14) 



which holds when [A, [A, B\\ = [B, [A, B]] = 0. Using Eqs. (JT3J) and (JUJ) one can verify the 
following relations: 



Tr(x[x, 


I x ,p\]) 


= o, 


Tr(V(x)[x, 


[ X ,P\]) 


= o, 


Tt(V 2 (x)[x, 


[ X ,P\]) 


= o, 


Tr(p x [x, 


[ x ,pW) 


= o, 


Tr (pl[x, 


i x iP]\) 


= -2h\ 


Tr(pl[x, 


[ x ,pW) 


= -l^ 2 (pl), 


V(x)p 2 x } [x, 


x ,p}}) 


= -Ah 2 (V(x)) 



(15) 

which are used to evaluate the effect of the environmental localization term. 

Given these expressions the time-evolution of the energy dispersion can be evaluated. 
Specifically: 

d(5E) 2 



dt 



d(H 2 ) 2(H) d(H) 

dt [ ' dt 
D (3(p 2 x )/m 2 + 2{V(x))/mj - 2 ([p 2 x )/2m + (V(x))) D/m 

2D(p 2 x )/m 2 , 



(16) 



where H is the system Hamiltonian given by H = p 2 ,/2m + V(x) and where we have used 
the result ^ that d(H)/dt = D/m. The extension to the two-dimensional H + H2 system 
is straightforward, giving 



d(5E) 2 2D 2 2D , 



(it ;U 2 31 /i 2 

and the rate of the average energy increase becomes 4D / \i 



-t(p 2 J + -t(p 2 v }, (17) 



V. DECOHERENCE IN ONE DIMENSIONAL SCATTERING 



A. Model Systems 

As is well known |25J] , systems with increasing nonlinearity are expected to show increasing 
quantum contributions to the dynamics. Hence we consider the effect of decoherence on 
several one-dimensional barrier potentials of different nonlinearities. Two different types of 
barriers are designed for each degree of nonlinearity, a single barrier potential (denoted SB) 
showing no resonance and a double barrier potential (denoted DB) that shows a typical 
resonance in the transmission probability versus initial energy. The insights gained from 
these simple systems help shed light on the H + H2 results. 

For the one-dimensional cases, the computational formulae introduced above reduce to 
the analogous equations for one degree of freedom. The numerical methods used to obtain 
transmission probabilities in the quantum and classical systems are the same as described 
in the H + H 2 reaction case, the only difference being that here the absorbing potential at 
the grid boundary is not implemented. Instead, a sufficiently large position range is used 
so that the wavepacket does not encounter the boundary during the time evolution. The 
labels QM, CM, QMD and CMD have the same meaning here as described in the H + H 2 
case. For the open quantum system a total of 500-1000 wavepackets are used, and for both 
the closed and open classical systems 10 4 trajectories are used. All variables below are in 
dimensionless units. 

The Hamiltonian for one- dimensional scattering is 
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with mass m and potential V(x). Here m = 1 and % = 0.1. The initial wavefunction is chosen 
as a minimum uncertainty Gaussian wave packet F(x) [Eq. (jHj)], and the associated Wigner 
function is used as the initial classical ensemble. In order to set SE of the initial wavepacket 
the same method as described above is used, except that here E v is zero, and the mass and 



HE 



h are different. We choose p x0 = (p x ) = \J2mE - h /47 2 , and thus 5E = J - 32 ^ m 2 + 2 ^ m - 
Specifically, we take 

Vsb(x) = V [cosh (ax)]" 2 , (19) 
V DB (x) = V SB (x -P) + V SB (x + (3) (20) 

as potentials for the SB case and for the DB case, respectively. The degree of system 
nonlinearity is estimated via a characteristic potential length, Xn{= \d x V /d' x l+1 V\ 1 / n ) [25J: 
the smaller the Xn, the larger the nonlinearity. Since the derivatives of these potentials obey 
the following relations: 

d 2n+1 V{x)/dx 2n+1 = (-Aa 2 ) n dV(x)/dx, (21) 
d 2n V(x)/dx 2n = (-Aa 2 ) n - 1 dV 2 (x)/dx 2 } (22) 

X2n+i = 1 /2a for both the SB and DB cases. Hence, increasing a leads to the larger 
nonlinearity. For the SB case, we examine two weakly nonlinear potential (denoted 

SBW) and a strongly nonlinear one (denoted SBS) with a = 0.5 and a = 10, respectively. 
For the DB case, (3 = 0.1 and only a strongly nonlinear case with a = 10 is examined. All 
Vq are chosen to giv e the same barrier height of 2.0 and numerical parameters are chosen to 
ensure accuracy [2fi| . 

Figure El shows the potential energy for SBW, SBS, and DB along with their first deriva- 
tives. Since the nonlinearity is inversely proportional to the width of the potential energy, a 
narrow/wide barrier has a large/small nonlinearity. This can be clearly seen in Fig. El Thus, 
noting that the first derivatives of SBS and DB are approximately 20 times larger than that 
of SBW, we anticipate that the quantum corrections for SBS and DB should be comparable 
to one another and larger than that of SBW and that the magnitude of the decoherence 
needed for the emergence of classicality in the SBW case should be smaller than that of the 
SBS case. This is indeed observed below. 
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B. Tunneling Regime 



Figure |U shows the results of the tunneling calculations. For cases with no decoherence, 
in addition to the classical (CM) and quantum (QM) results, we show the exact quantum- 
mechanical transmission probabilities for the single barrier potentials, denoted QMA. They 
are given by [27] 



Pr = smh 2 (7rp x0 / ha) /[smh. 2 (7rp x0 /ha) + cos1i 2 {7t/2a/c 2 — 1}] (23) 

where c 2 = 8mVo/h 2 a 2 > 1, and used as a computational check. In the absence of 
decoherence the strongly nonlinear cases, SBS and DB, show larger differences between the 
quantum and classical results than does the weaker nonlinearity case. Further, tunneling in 
the quantum cases is clear, as is the broad resonance around E = 2.6 in the double barrier 
case, reflecting the presence of a metastable state in the quantum double barrier. 

In Fig. 0J the decoherence shown is of magnitude that lead to roughly good agreement 
between quantum and classical results. In panel (a) this corresponds to D = 6 x 10~ 4 , 
but in panels (b) and (c) the required D = 3 x 1CT 2 , i.e. it is 50 times larger than that 
in (a). A comparison of results between the closed and open systems shows that introduc- 
ing decoherence substantially improves quantum-classical correspondence. Several relevant 
observations are in order. First, the quantum resonance in the DB case is suppressed by 
increasing decoherence, leading to better agreement between quantum and classical trans- 
mission probabilities. This is discussed further in Section IV CI Second, in order to reach 
classicality, the decoherence needed for both the SBS and DB cases is 50 times larger than 
that of the SBW case, due to the higher nonlinearity of their potentials. It is tempting to 
envision that classicality in the double barrier case might be reached with smaller decoher- 
ence than that of the single barrier SBS case of the same nonlinearity, since the quantum 
DB results shows a resonance that is absent in the SBS case. This view comes from the 



semiclassical picture |28[ suggesting that the resonance region may be more sensitive to 
decoherence than is the threshold region since resonance require constructive interferences 
between trajectories. However, careful examination shows that it is the nonlinearity that 
dictates the magnitude of the decoherence required to reach classicality. 

The third observation requires more attention. Note that Fig. 0] captures an overall 
qualitative picture of decoherence effects on the transmission probabilities. That is, in both 
the quantum and classical cases, decoherence enhances the transmission at energies below 
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the barrier height and suppresses it above the barrier height, similar to the trend seen in 
H + H2. This behavior can be understood by considering the increase in energy width 8E 
associated with introducing decoherence of the type given in Eq. JI}. Figure shows the 
computed growth of 5E for the closed and open systems in the SBS case. The quantum and 
classical energy widths are seen to be in good agreement as a function of time, and become 
comparable to the 2-unit high barrier. Noting that the average energy (E) for the open 



121 as well, one can understand from the 



system increases linearly with the rate of D/m 
classical point of view why decoherence enhances the transmission at E = 1.4, suppresses 
it at E = 3.0, and has little effect at E = 2.0 [Fig. E]^b)]. Specifically, at the collision time 
of t ~ 3, in the presence of decoherence, SE is approximately 1.05, 0.75, and 0.70 [Fig. |S] 
for the initial average energy E— 3.0, 2.0, and 1.4, respectively, and is much larger than the 
average energy increase (~ 0.06). Consequently, when the wavepacket of E = 3.0 strikes the 
barrier, the decoherence-induced SE increases that part of the wavepacket below the barrier 
height, increasing the reflected component, and reducing the transmission probability. (Note 
that although the average energy increases, decoherence-induced changes in the portion of 
the wavepacket above the barrier will not alter the transmission, since it will transmit in any 
case.) By contrast, for the case of E = 1.4, decoherence increases the portion of the energy 
above the barrier, allowing it to pass over the barrier, thus increasing the transmission 
probability. For the case of E = 3.0 and E = 1.4 the same trends are seen in both the 
quantum and classical results. On the other hand, for E = 2, equal to the barrier height, 
the transmission probability is rather insensitive to the decoherence. Similar observations 
apply to the other barrier cases. 



C. Resonance Regime 

Consider now the resonance region. This region has a characteristic width denoted 5Er. 
It is reasonable to assume that when D is sufficiently large to ensure that 5E is on the order 
of 5Er at the time of the scattering resonance then the resonance feature is eliminated. 
To obtain a crude estimate of the D for which this occurs consider Eq. (|T6*|) : 

^£ = 2D (p 2 x )/m 2 « ADEo/m, (24) 
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where Eq is (H) at time zero. Here, in order to obtain a simple estimate, the initial kinetic 
energy has been assumed constant in time. Noting that here SE at time zero is ignorable 
compared to the growth rate, gives 



A comparison with Fig. El shows that this drastic approximation simply replaces the satu- 
ration behavior by linear growth of 5E, overestimating the growth of 5E. 

The simplest D estimate is then that the resonance width SEr, which occurs at time £r, 
will be washed out for decoherence strengths Dr on the order of: 



Sample computations using this estimate are shown in Fig. |U] for the DB case with 
varying mass m (which implicitly changes SEr and In the cases shown, SEr varied 
from 0.6 to 0.3 and Dr from 0.007 to 0.0035. In all cases we compare the decoherence free 
result (denoted QM) to the result with either decoherence with Dr or, for comparison, with 
two different values of D, one being Dr. As is evident from Figure |B1 the characteristic 
resonance shape disappears with the application of decoherence of magnitude Dr in each 
case. Interestingly, as in H + H 2 , the resonance region disappears via a lowering of the initial 
peak that indicates resonant behavior, leaving a smooth falloff of Pr with decreasing energy. 
The smaller values of D, shown in panels (a) and (b) for comparison, are insufficiently large 
to completely eliminate the resonance. 

D. H + H 2 

This argument can be extended to the resonance region in H + H2. Since our quantum 
H + H2 calculation does not provide the increase of SE (due to absorbing boundary con- 
ditions), SE is estimated from the classical ensemble. At t = 150 fs ~ tf/2 and at E = 
0.9 eV, which is near the resonance dip, the smaller decoherence magnitude used above 
(D = 2.47 x 10~ 35 kg-J/s) gives an average energy increase of ~ 0.004 eV, which is relatively 
small compared to the average energy, and an energy width SE ~ 0.08 eV, which is roughly 
comparable to the resonance width. Hence, one would expect only partial elimination of 
the resonance, which is indeed the case [Fig. Q]. By comparison, the larger value of D used 




(25) 



D R « (SER) 2 m/[4E t R ] 



(26) 
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[Fig. |2] is expected to have a 5E that is roughly 0.4 eV and an average energy increase of 
~ 0.09 eV, giving the extremely broad featureless results in Fig. El The resultant behavior 
in the H + H 2 case, is then entirely consistent with the broadening due to the decoherence 
induced SE. 



VI. SUMMARY 



Decoherence has been shown to significantly alter both the quantum and classical H + H 2 
reaction probabilities, leading to quantum -classical agreement at realistic decoherence val- 
ues. Unexpected behavior in the tunneling region, where the reaction probability is increased 
below the reaction threshold, but decreased above the threshold, lead to a more detailed 
examination of decoherence in the tunneling regime of simple one dimensional problems. 
The results of this study showed that (a) the extent of potential nonlinearity dictated the 
extent of decoherence required to reach quantum-classical agreement: that is, the larger the 
nonlinearity, the higher the required decoherence parameter, and (b) that the increase in 
energy width, due to environmental-induced localization, was responsible for the observed 
correspondence in probabilities in the tunneling and resonance regimes. 

Studies of this type also made clear the general dearth of information on actual decoher- 
ence parameters for realistic systems. For this reason, realistic studies of this type are in 
progress 
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Appendix: Values of D 



Within Caldeira-Leggett type models, the constant D is related to the relaxation constant, 
/by |3o| 

D = 2fmkT 1 (27) 

where m is the system mass, k the Boltzmann constant, T the temperature of the envi- 
ronment, and / the strength of the system-environment coupling. Since a typical velocity 
relaxation time for small molecules is 0.01 - 0.1 ps [31], D typica i is roughly 10~ 33 — 10~ 34 
kg-J/s for the mass of hydrogen at room temperature. This gives, for our choice of the 
smaller value of D that D/D ty p ica i ~ 10~ 2 — 10 _1 . 
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FIGURE CAPTIONS 



Figure 1. Collinear reactive transition probabilities in H + H2 versus the initial total energy 
E, without and with decoherence. Here D = 2.47 x 1CT 35 kg-J/s. 

Figure 2. As in Fig^ but with 20 times stronger decoherence, i.e. , D = 4.94 x 10~ 34 kg-J/s. 

Figure 3. Left panel: Potential energy of SBW with scale at left, and the first derivative of 
SBW with scale at right. For comparison, the much sharper SBS potential energy is shown 
as well. Right panel: Potential energies and their first derivative for the SBS and DB cases. 
Thin lines and thick lines correspond to SBS and DB, respectively. 

Figure 4. Transmission probabilities as a function of initial energy, E, without and with 
decoherence. (a) SBW case, (b) SBS case, and (c) DB case. All variables are in dimensionless 
units. 

Figure 5. Time dependence of 5E for the open and closed SBS system. Solid lines and cross 
symbols denote the quantum and classical results, respectively, for the open system; from 
top to bottom E takes the values 3.0, 2.0 and 1.4. A thick dashed line along SE= 0.1 denotes 
the quantum and classical results for the closed system. All variables are in dimensionless 
units, and all initial SE ps 0.1. 

Figure 6. Transmission probabilities for the DB case, with and without decoherence, for 
various masses, to test the predicted Dr [Eq. (J2HJ)]- (a) m = 1, £r=6, 5Er=0.6, Dr = 
£2=0.006, D!=0.003, (b) m = 2, t R =8, 6E R =0.5, D R = D 2 =0.007, £>i=0.002, (c) m = 3, 
t R =9, 5E R =0.3, D R = £>=0.0035. 
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